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We introduce a framework of opto-mechanical systems that are driven with a mildly amplitude-modulated 
light field, but that are not subject to classical feedback or squeezed input light. We find that in such a system 
one can achieve large degrees of squeezing of a mechanical micromirror - signifying quantum properties of 
opto-mechanical systems - without the need of any feedback and control, and within parameters reasonable 
in experimental settings. Entanglement dynamics is shown of states following classical quasi-periodic orbits 
in their first moments. We discuss the complex time-dependence of the modes of a cavity-light field and a 
mechanical mode in phase space. Such settings give rise to certifiable quantum properties within experimental 
conditions feasible with present technology. 



Periodically driven quantum systems exhibit a rich behavior 
and display non-equilibrium properties that are absent in their 
static counterparts. By appropriately exploiting time-periodic 
driving, strongly correlated Bose-Hubbard-type models can 
be dynamically driven to quantum phase transitions (TJ, sys- 
tems can be dynamically decoupled from their environments 
to avoid decoherence in quantum information science 0, and 
quite intriguing dynamics of Rydberg atoms strongly driven 
by microwaves [3| can arise. It has also been muted that 
such time-dependent settings may give rise to entanglement 
dynamics in oscillating molecules J4) . A framework of such 
periodically driven systems is provided by the theory of lin- 
ear differential equations with periodic coefficients or inho- 
mogeneities, including Floquet's theorem Q. 

In this work, we aim at transferring such ideas to describe 
a new and in fact quite simple regime of opto-mechanical 
systems, of micromirrors as part of a Fabry-Perot cavity 
E] 13 El): So to one of the settings iflOl ITT1 [Til fl3l fl4l 
that are the most promising candidates in the race of explor- 
ing certifiable quantum effects involving macroscopic me- 
chanical modes. This is an instance of a regime of driving 
with mildly amplitude-modulated light. We find that in this 
regime, high degrees of squeezing below the vacuum noise 
level can be reached, signifying genuine quantum dynam- 
ics. More specifically, in contrast to earlier descriptions of 
opto-mechanical systems with a periodic time-dependence in 
some aspect of the description, we will not rely on classical 
feedback based on processing of measurement-outcomes - a 
promising idea in its own right in a continuous-measurement 
perspective {15\ [l6l - or resort to driving with squeezed light. 
Instead, we will consider the plain setting of a time-periodic 
amplitude modulation of an input light. The picture developed 
here, based in the theory of differential equations, gives rise to 
a framework of describing such situations. We find that large 
degrees of squeezing can be reached (complementing other 
very recent non-periodic approaches based on cavity-assisted 
squeezing using an additional squeezed light beam ifTTl ). It 
is the practical appeal of this work that such quantum signa- 
tures can be reached without the necessity of any feedback, 
no driving with additional fields, and no squeezed light input 
(the scheme by far outperforms direct driving with a single 
squeezed light mode): In a nutshell, one has to simply gen- 



tly shake the system in time with the right frequency to have 
the mechanical and optical modes rotate appropriately around 
each other, reminding of parametric amplification, and to so 
directly certify quantum features of such a system. 

Time-dependent picture of system. Before we discuss the 
actual time-dependence of the driven system, setting the stage, 
we start our description with the familiar Hamiltonian of a 
system of a Fabry-Perot cavity of length L and finesse F being 
formed on one end by a moving micromirror, 

H = huu c a) a + -hui m (p 2 + q 2 ) — fiG^a) aq 

oo 
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Here, ui m is the frequency of the mechanical mode with 
quadratures q and p satisfying the usual commutation rela- 
tions of canonical coordinates, while the bosonic operators a 
and are associated to the cavity mode with frequency co c 
and decay rate k = nc/(2FL). Gq — \J hj (muo m ) uj c / L is 
the coupling coefficient of the radiation pressure, where m is 
the effective mass of the mode of the mirror being used. Im- 
portantly, we allow for any periodically modulated driving at 
this point, which can be expressed in such a Fourier series, 
where Q, — 2it/t and t > is the modulation period. The 
main frequency of the driving field is ujq while the modulation 
coefficients {E n } are related to the power of the associated 
sidebands {P n } by \E n \ 2 = 2nP n /(hu)Q). The resulting dy- 
namics under this Hamiltonian together with an unavoidable 
coupling of the mechanical mode to a thermal reservoir and 
cavity losses gives rise to the quantum Langevin equation in 
the reference frame rotating with frequency ui , q = uj m p, and 

p = -uj m q-j m p+G a^a + ^, (2) 

OO 

a = -(«+ iA )a + iG a aq+ ^ E n e~ innt + V2~K,a in , 
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where Ao = uj c — ujq is the cavity detuning. j m is here an 
effective damping rate related to the oscillator quality factor 
Q by 7, Tl = u> m /Q. The mechanical (£) and the optical (a m ) 
noise operators have zero mean values and are characterized 
by their auto correlation functions which, in the Markovian 
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approximation, are 

<£(W)+£(*')£(*)>/2 = 7m(2n+ l)6(t-t') (3) 

and (a m (i)a in t(t')) = 5(t-t'), where n = [cxp(|^f 
is the mean thermal phonon number. Here, we have implic- 
itly assumed that such an effective damping model holds ifTHl . 
which is a reasonable assumption in a wide range of parame- 
ters including the current experimental regime. 

Semiclassical phase space orbits. Our strategy of a solu- 
tion will be as follows: we will first investigate the classical 
phase space orbits of the first moments of quadratures. We 
then consider the quantum fluctuations around the asymptotic 
quasi-periodic orbits, by implementing the usual linearization 
of the Heisenberg equations of motion ifTTl fT2l (excluding the 
very weak driving regime). Exploiting results from the theory 
of linear differential equations with periodic coefficients, we 
can then proceed to describe the dynamics of fluctuations and 
find an analytical solution for the second moments. 

If we average the Langevin equations Q, assuming 
(a^a) ~ |(a)| 2 , (aq) ~ (a)(q) (true in the semi-classical driv- 
ing regime we are interested in), we have a nonlinear differen- 
tial equation for the first moments. Far away from instabilities 
and multi-stabilities, a power series ansatz in the coupling Go 
(0}(t) — Y^nLoOjifyGo' is justified, where O = a,p,q. 
If we substitute this expression into the averaged Langevin 
equation we get a set of recursive differential equation 
for the variables Oj(.). The only two nonlinear terms in Eq. 
Q are both proportional to Go, therefore, for each j, the dif- 
ferential equation for the set of unknown variables Oj(.) is a 
linear inhomogeneous system with constant coefficients and 
r-periodic driving. Then, after an exponentially decaying ini- 
tial transient (of the order of l/7 m ), the asymptotic solutions 
Oj will have the same periodicity of the modulation [5 1, jus- 
tifying the Fourier expansion 
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Substituting this in Eq. d2), we find the following recursive 
formulae for the time independent coefficients O n j, 



q n ,o = Pn.a = 0, 

a n . = £7_ n /(« + i(Ao + nfi)), 



(5) 
(6) 



corresponding to the zero coupling limit Go = 0, while for 
each j > 1, we have 
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Within the typical parameter space, considering only the first 
terms in the double expansion (HJ, corresponding to the first 



side bands, leads to a good analytical approximation of the 
classical periodic orbits, see Fig. 1. On physical grounds, this 
is expected to be a good approximation, since Go <C w m , 
and because high side-bands (of frequency n£l) fall outside 
the cavity bandwidth, nil > 2k. What is more, the decay 
behavior of E n related to the smoothness of the drive inherits 
a good approximation in terms of few sidebands. 
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FIG. 1: Phase space trajectories of the first moments of the mirror (a) 
and light (b) modes. Numerical simulations for t € [0, 50r] (black) 
and analytical approximations of the asymptotic orbits (green). 

Quantum fluctuations around the classical orbits. We will 
now turn to the actual quantum dynamics taking first mo- 
ments into account separately when writing any operator as 
0(t) = (0)(t) + 50(t). The frame will hence be provided 
by the motion of the first moments. In this reference frame, 
as long as |(a)| ^> 1, the usual linearization approximation to 
Q can be implemented. In what follows, we will also use the 
quadratures 5x = (8a+5aJ)/V2and5y — —i(5a—8a))/*j2, 
and the analogous input noise quadratures x m and y m . For the 
vector of all quadratures we will write u = (5q, Sp, Sx, 5y) T , 
with n = (0, £, y/2KX in , ^/2ny m ) T being the noise vector 
IfTTl [18 1 . Then the time-dependent inhomogeneous equations 
of motion arise as u(t) = A(t)u(t) + n(t), with 



A(t) 
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G x (t) G y (t) 

-k A(t) 
-A(t) -K 



(8) 



where the real A(t) contains the time-modulated coupling 
constants and the detuning as G(t) = G x (t) + iG y [t), 



G(t) = V2(a(t))G , A(t) = A - G (q(t)). 



(9) 



From now on we will consider quasi-periodic orbits only, so 
the long-time dynamics following the initial one, when the 
first moments follow a motion that is r-periodic. Then, A is 
r-periodic, and hence 



A(t)=A(t + r) 



J2 A n e mnt . 

71 — - — OO 



(10) 



In turn, if all eigenvalues of A(.) having negative real parts 
for all t £ [0,r] is a sufficient condition for stability. From 
the Markovian assumption Q, we have 



{ni(t)nj{t') + nj(t')n t (t))/2 = 6(t - t') D % 



(11) 
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where D = diag(0, 7 m (2n + 1), k, k). The formal solution of 
Eq. © is El 

u(t) = U(t,t )u(t ) + [ U(t,s)n(s)ds, (12) 

where U(t,to) is the principal matrix solution of the ho- 
mogeneous system satisfying U(t,to) = A(t)U(t,to) and 
C/(to, to) — 1. From Eqs. ([8] [12}, we have as an equation 
of motion of the covariance matrix (CM) of the two modes 

V(t) = A(t)V(t) + V(t)A T (t) + D. (13) 



Here, the CM V(.) is the 4 x 4-matrix with components 
Vi.j = (uiUj + UjUi) /2, collecting the second moments of 
the quadratures. This is again an inhomogeneous differential 
equation for the second moments which can readily be solved 
using quadrature methods, providing numerical solutions that 
will be used to test and justify analytical approximate results 
in important regimes. Moreover, now the coefficients and not 
the inhomogeneity are r-periodic, A(t) = A(t + r). Again, 
we can invoke results from the theory of linear differential 
equations to Eq. ( p"3j ) [5I|: We find that in the long time limit, 
the CM is periodic and can be written as 



V{t) = Y,V n e 



i nil! 



(14) 



An analytical solution for V(.), is most convenient in the 
Fourier domain, f(u>) = J^°° e~ luJt f(t)dt, giving rise to 

oo 

— iuiu(u>) + A„{t(cj — nfl) = — n(u)) . (15) 



If A n9 LQ — 0, corresponding to no-modulation, we are in 
the usual regime where the spectra are centered around ±u> m 
for the mechanical oscillator and around ±A for the optical 
mode. The modulation introduces sidebands shifted by ±nf2. 
If the modulation is weak, only the first two sidebands at ±f2 
significantly contribute. For strong modulation also further 
sidebands play a role: Disregarding higher sidebands means 
truncating the summation to ±JV, (valid if u(ui ± Nil) ~ 0). 
Then Eq. ( 15 i can be written as A{uj)u(lo) = n(u>), where 



u T (to) = (Wiuj - Nn), . . . ,u t (lo), . . . ,u I (u + NQ)) and 
n T (cu) = (n T '(u - Nil), . . . , n T '(«), . . . , n T '(w + NQ)) are 
4 x (2A^ + 1) vectors, while, in terms of 4 x 4 blocks, 



A{w) 



A-i 



A. 



A- 



A\ B_(jv-i) A-\ 

A2 A\ _B_(]v_2) 



A2N 

with B k = A - i{u + Ml). 



A! Br, 



(16) 



From the correlation properties of the noise vector n(.), we 
have that 



<pij(u),oj') = (ni(w)n^(w') +rij(uj')ni(oj))/2 



2N 



= 8{u) — {J — n£l)D n , 



(17) 
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where Do = diag(£>, D, . . . , D), then D\ is the matrix that 
has D on all first right off diagonal blocks, D2 on the sec- 
ond off diagonals, with D n analogously defined, and -D_ n = 
D^- We now define the two frequency correlation func- 
tion as Vi,j{w, w') = (ui(w)u*(u>') + u*(oj')ui(uj))/2. We 
have V{uj,J) = A _x (o;)</»(a;,a; / )[A- 1 (a;')]t. We are in- 
terested only on the central 4x4 block of V which we call 
V(ui, u>') = [V(oj, u>')]4. From <p(ui, u>'), we find 



2N 
n=-2N 



nfl), (18) 



where V n {uo) = [A~ 1 (ui)D n [A^ 1 (uj — nCl)y]^. This means 
that the driving modulation correlates different frequencies, 
but only if they are separated by a multiple of the modulation 
frequency fl By inverse Fourier transforms we recover the 
time periodic expression for the CM, where the all the com- 
ponents V n are given by the integral of their noise spectra, i.e., 



V n = 



1 

27 



+OO 



V n {uj)duj. 



(19) 



Squeezing and entanglement modulation. We will now 
see that the mild amplitude-modulated driving in the cool- 
ing regime is exactly the tool that we need in order to ar- 
rive at strong degrees of squeezing, in the absence of feed- 
back or squeezed light. We will apply the previous general 
theory to setting of an optomechanical system that is experi- 
mentally feasible with present technology. In fact, all values 
that we assume have been achieved already and reported on in 
publications with the exception of assuming a relatively good 
mechanical Q-factor. The reasonable set of experimental pa- 
rameters [9| that we assume is L = 25mm, F = 1.4 x 10 4 , 
uj m = 2ttMHz, Q = 10 6 , m = 150ng, T = 0.1K. We 
then consider the - in the meantime well known - self-cooling 
regime [7] in which a cavity eigen-mode is driven with a red 
detuned laser A — tu m (with wavelength A = 1064 nm), but 
we also add a small sinusoidal modulation to the input am- 
plitude with a frequency il = 2ui m , so twice the mechanical 
frequency. To be more precise we choose the power of the 
carrier component equal to Po = lOmW, and the power of the 
two modulation sidebands equal to P±\ = 2m W. 

We approximate the asymptotic classical mean values in 
Eq. Q by truncating the series only to the first terms with 
indexes j = 0, . . . ,3 and n = — 1, 0, 1, giving rise to a good 
approximation. Fig. [T] shows that, after less than 50 modu- 
lation periods, the first moments reach quasi-periodic orbits 
which are well approximated by our analytical results. 
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FIG. 2: (a) Variance of the mirror position and (b) light-mirror 
entanglement En as functions of time. In both (a) and (b) the 
non-modulated driving regime (blue), the modulated driving regime 
(green) and the numerical solutions (black dashed/dotted) are plotted, 
(a) also shows the standard quantum limit (red dashed) at 1/2, the 
minimum eigenvalue of the mirror covariance matrix (black dashed) 
and its analytical estimation p,0\ in the RWA (orange). 



Conclusions and outlook. In this work we have introduced 
a framework of describing periodically amplitude-modulated 
optomechanical systems. Interestingly, such a surpringly sim- 
ple setting feasible with present technology |9l leads to a set- 
ting showing high degrees of mechanical squeezing, with no 
feedback or additional fields needed. We hope that such ideas 
contribute to experimental studies finally certifying first quan- 
tum mechanical effects in macroscopic mechanical systems, 
constituting quite an intriguing perspective. 
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In order to calculate the variances of the quantum fluctua- 
tions around the classical orbits, we truncate the sum in Eq. 
( 1 5 i to N — 2 and we apply all the previous theory to find the 
covariance matrix V. In Fig. [2] we compare two regimes: with 
or without {P±\ — 0) modulation (computed analytically and 
numerically). We see that the modulation of the driving field 
causes the emergence of significant true quantum squeezing 
below the Heisenberg limit of the mechanical oscillator state 
and also the interesting phenomenon of light-mirror entan- 
glement oscillations. This dynamics reminds of the effect of 
parametric amplification 01311161 . as if the spring constant of 
the mechanical motion was varied in time with just twice the 
frequency of the mechanical motion, leading to the squeezing 
of the mechanical mode. For related ideas of reservior engi- 
neering, making use of bichromatic microwave coupling to a 
charge qubit of nano-mechanical oscillators, see Refs. Il23l . 
Here, it is a more complicated joint dynamics of the cavity 
field and the mechanical mode - where the dynamics of the 
first and the second moments can be separated - which for 
large times yet yields a similar effect. Indeed, this squeezing 
can directly be measured when considering the output power 
spectrum, following Ref. [20|, and no additional laser light is 
needed for the readout, giving hence rise to a relatively sim- 
ple certification of the squeezing. Entanglement here refer 
to genuine quantum correlations between the mirror and the 
field mode, as being quantified by the logarithmic negativity 
defined as .E/v(p) — log ||p r ||i, essentially the trace-norm of 
the partial transpose, which is a proper entanglement measure 
ET1 |22"1 . The minimum eigenvalue of the mirror covariance 
matrix - the logarithm thereof typically referred to as single 
mode squeezing parameter - is almost constant and this means 
that the state is always squeezed but that the squeezing direc- 
tion continuously rotates in phase space with the same period 
of the modulation. Calling this rotating squeezed quadrature 
5x ij, a rough estimate of its variance can be calculated in the 
rotating-wave approximation (RWA, compare, e.g., Ref. G4l0 . 

2 1 2 K (Go-G_ 1 )(G n + G_ 1 (n+l)) /om 

{6XR} ~2 +n ~ (7 m + 2«)(Gg-G?. 1 + 2 7m «) (20) 
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APPENDICES 

In this appendix, we will summarize a number of additional 
statements useful for the main text. 



Differential equations 

In this appendix, we will consider several issues related to 
solutions of inhomogeneous first-order differential equations 
that are used in the main text. Most of the presented material is 
standard material from the theory of time-dependent periodic 
linear systems. At the end, yet, we will consider a bound to the 
error made when assuming a periodic solution that is tailored 
to the specific situation studied above. Consider a linear first- 
order system, 



±(t) = B{t)x{t), 



(21) 



where B(t) being some complex square matrix with entries 
dependent on t > to with some to > 0. Then the linear first- 
order system has a unique solution for 2; (to) = £0 f° r a ll times 
t > to- The principal matrix solution is the solution of 



P(t,t ) = B(t)P(t,t ), 



(22) 



with P(t,to) = 1. Floquet's theorem now states the follow- 
ing: 

Lemma 1 (Floquet's theorem) If B(.) is periodic, B(t) = 
B(t + r) for some r > for all t > to, then the principal 
matrix solution has the form 



P(t,to) = X(t,to)e [t - to)Y{to \ 



(23) 



where the matrices X(., .) and Y(.) are T-periodic in all their 
arguments and X (to, to) = 1. 

So if one factors out an exponential, the remainder is periodic 
in time. The eigenvalues of the monodromy matrix 



M (t ) = P(t + t, to) 



(24) 



are then referred to as Floquet multipliers, and the eigenvalues 
of Y(to) are known as Floquet exponents. Floquet exponents 
are t independent, in fact the matrix Y(t ) is similar to Y (t' ) 
for any t' Q . If all Floquet exponents have a negative real part, 
the system is asymptotically stable. 



Lemma 2 (Solution to inhomogeneous problem) The solu- 
tion to the inhomogeneous system 

x(t) = B(t)x(t)+g(t), (25) 

with initial condition a; (to) = %o is given by 

x(t) =P(t,t )x + f dsP(t,s)g(s). (26) 

Jto 

We will now show that under simple conditions that, when 
both B(.) and <?(.) are T-periodic, we will asymptotically ar- 
rive at a solution with the same time period. Let us define the 
interval I = [0, r]. In order to prepare that statement, we will 
use the following bound. 

Lemma 3 (Bound from Floquet exponents) In the above 
notation, if the system is stable and ift — to > 1, tlien 

max||e( t -* 0+tl ) Y ( tQ - ll )|| < cn{t-to+T) n - 1 e x{ - t - t °\ (27) 



where 



max||W(u)|| ||W - V)II> 



(28) 



where W(u) is a similarity transformation that brings Y{u) 
to a Jordan normal form. 



Here, 



A = max re(Aj), 



where Xj are the Floquet exponents. The norm 



||A|| = sup 



Nl 



(29) 
is the norm 
(30) 



induced by the usual Euclidean vector norm ||.||. So up to 
a constant, the convergence is essentially governed by the 
largest real part of the Floquet exponents. To show this bound, 
note that 

e (t-t Q +u)Y(t -u) = w -\to - u) [® 1 M J (5)e^ s ] W(t - u), 

(31) 



where 



Mj(6) = 



1 5 g ••• 
1 5 f 



K-i)! 



S = t-to+u, (32) 



and rij is the dimension of the Jordan block associated with 
the Floquet exponent Xj, giving rise to the above expression 
( 27 1, by acknowledging that 



\\Mj(5)\\ < rijS^- 1 < n5 n 



(33) 



for all j, by bounding the ||.|| norm from above by the ||.||2- 
norm, and since S > 1. We can now turn back to our original 
problem of the periodicity of the asymptotic solution. 
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Theorem 1 (Asymptotic periodicity) If tlie system is stable 
and both B{.) and <?(.) are T-periodic (including the case 
of constant B or g), then the solution x(.) of the inhomoge- 
neous problem is asymptotically also T-periodic, with known 
bounds, in time exponentially suppressed, on the error made. 



From ( 26 1, with the new integration variable u = s — r, one 



can write for the solution of the inhomogeneous system 

x(t + r) = P(t+T,t )x (34) 



duP(t + r, u + r)g(u + r) 



to — T 



P(t + r, t )x 



duP(t, u)g(u), 



to—T 



using the appropriate periodicities of g and P. Now we try 
to bound the following quantity, which should exponentially 
decay to zero in order to prove the theorem, 



i{t + r)-x{t) = X{tMW Y(ta) + l]e {t ' ta)Y[to) x 



duP(t, u)g(u). 



t a — T 



With a new variable v — to — u, the integral in Eq. ( [35 
becomes 



dvX(t,t Q - v)e {t - t0+v)Y{to ~ v) g(t - v), 



which is bounded from above by 



rm max \\g(v) || max || e (*-*°+«)n*o-«) | 



where 



m= max \\X(t,t')\\ 

(t,t')<£lxl 



(36) 



(37) 



(38) 



Now, by using the triangular inequality, we find the following 
bound for the norm of the vector in i 



\\x(t + t) - x(t)\\ < mmax||e (t - to+1,)Y(to - t ' ) || 



2\\x \\ + Tmax||#(» 



(39) 



Finally, by using Lemma [3] we have 

\\x(t + t) - x(t)\\ < e^-^mcnit-to + T) 71 - 1 

x(2\\x \\ +Tmax\\g(v)\\). (40) 

\ V&I / 

Now it is clear that, in the long time limit (t — t — > oo), the 
first factor e A ( t_t °) exponentially suppresses the RHS of Eq. 



Squeezing in the rotating wave approximation 

Suppose that we are in the particular resonance condition 
such that 



A(t) 



0J„ 



A n e mQt 



(41) 



and = 2u> m . Then, in the limit of cj m ^> |G(<)|, K, a sim- 
ple analytical expression for the covariance matrix can be ob- 
tained in the rotating wave approximation. We observe that we 
have two arbitrary degrees of freedom that we can choose to 
simplify the calculation: The global phase of the driving laser 
and the initial time corresponding to t = 0. If we expand the 
coupling coefficient G(t) = X)^L-co G n e mflt , then we can 
assume without loss of generality assume that Go, G_i € ML 
Now we move to an interaction picture introducing the 
slowly varying bosonic operators, 



a s e 



\b s = be 1 ' 



(42) 



In this reference frame, if we neglect terms rotating at fre- 



quency 2ui m , we obtain an equation analogous to Eq. ( 13 1 



where, 



A. 



and 



V s {t)=A s V s + V s A T s +D s (43) 



"7m G_!-G 

"7m G_!+G 

G_i - G -2k 0) 
G_i + Go -2k 

(44) 



D = diag( 7 (n + 1/2), 7 (n + 1/2), k, k). 



(45) 



We observe that, in the RWA, only the coefficients Go and 
Gi matter. They correspond to the cooling and the heating 
sidebands of the input driving laser. 



The stability condition for the differential Eq. (43 i is G_ x 



Gq < 2 7to k, which is always satisfied if the cooling process 
is predominant with respect to the heating. Differently from 
Eqs. ( [13) , Eq. ( |43] l has constant coefficients and therefore, in 
the stable regime, V s (t) reaches the an asymptotic constant 
value V s = lim^oo V s (t). 

The matrix V s can be calculated by imposing the derivative 



( 40 1 proving the theorem. 



in (43 1 equal to zero and solving the remaining linear system. 
We report only the mirror variances: {y s )\,x — /-, (V^)a,2 = 
/+ and (V s )\ : 2 = 0, where 

1_ 2^(Go±G_ 1 )[G nTG_ 1 (n+l)] 
J± 2 +n ( 7m + 2 K )(G2-G 2 _ 1+ 2 7mK ) ■ {0> 

Three particular limits. If G_i = 0, which corresponds to 
a red detuned driving laser without modulation, we observe 
the usual cooling of the mirror. For G_i = Go, we recover 
the QND measurement setting studied in Ref. lfT5ll . where a 
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symmetric driving with opposite detuning couples the light 
only with a quadrature of the mirror. Indeed, the variance V s ii 
is unaffected by the driving while the other is increased due to 
radiation pressure noise. Finally if G_i < Go, we observe 
that the mirror can reach a steady state which is both cooled 
and squeezed without the need of any feedback and control. 



Squeezed environments versus modulation 

In recent work [17] it has been shown that it is possible 
to squeeze the mechanical mode by using a constant driving 
field with A = io m and an additional squeezed vacuum beam 
with frequency resonant with cavity uj s — lo c . In the RWA 
formulation used in the previous section, this corresponds to a 
coefficient matrix A' having the same structure as in Eq. (44 1 
but with G = G' and G_i =0, while 



D' = diag( 7 (n + 1/2), 7 (n + 1/2), an, s _1 k), (47) 

where s is a parameter quantifying the squeezing of the en- 
vironment. If we solve the system for the mirror variances, 



instead of Eq. (46 1 we get 



f± 



K G 2 (2n- 



( 7m + 2 K )(G 2 + 2 lm n)' 



(48) 



Now we show that, within the validity of the RWA, for what 
concerns the squeezing of the mirror, this approach is equiv- 
alent to our scheme in which the cavity is driven by a single 
amplitude modulated laser and no squeezing is required. In 
fact, by choosing the modulation parameters such that 



Go = Gis- 1 ' 2 



,1/2 
,1/2 



)A 
)A 



(49) 
(50) 



the previous formula given in (46 1 reduces exactly to Eq. (48 1. 
This means that the two approaches are formally equivalent as 
far as the squeezing is concerned and the choice between one 
or the other depends only on technical factors. 

Finally, we may ask if a good mechanical squeezing is pos- 
sible by just driving the cavity with a single squeezed non- 
modulated field with A = uj m . Differently from the proposal 
of Ref. ifTTl . in this case the squeezed input noise is not res- 
onant with the cavity mode but rotates with a frequency ui m . 
This means that, in the RWA, the squeezing is averaged to 
zero and the optical environment looks like an effective ther- 
mal bath. Therefore, the only convenient choice is to use an 
additional squeezed beam which is resonant with cavity as in 
Ref. IfTTl . Only in this way the squeezing direction of the en- 
vironment does not rotate and can have a significant effect on 
the system dynamics. 



